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Based on the developed quantum microscopic theory, the interaction of weak electromagnetic 
radiation with dense ultracold atomic clouds is described in detail. The differential and total coop- 
' erative scattering cross sections are calculated for monochromatic radiation as particular examples 

^SJ , of application of the general theory. The angular, spectral, and polarization properties of scattered 

light are determined. The dependence of these quantities on the sample size and concentration of 
atoms is studied and the influence of collective effects is analyzed. 
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I. INTRODUCTION 



(~| I A great number of experimental and theoretical papers devoted to the study of cold and ultracold atomic 

. ensembles are related to the unique physical properties inherent in these systems and their broad possible 

^ , ' practical applications. Improvements in techniques for cooling atomic gases in atomic traps ^3 make their 

, ', use very promising in metrology, frequency standards, and quantum information problems Most of the 

methods for diagnostics of the states of these ensembles are based on their interaction with electromagnetic 
radiation. More importantly, the largest number of the many applications envisioned for ultracold atom 
^ ■ ensembles have at their foundation such interactions. 

\ The interaction of radiation with cold atomic clouds reveals a number of specific features. First, atoms cooled 

down to low temperatures have very high resonance scattering cross sections, and therefore multiple scattering 
should be taken into account even in quite rarefied clouds. Second, atomic clouds in traps have nonequilibrium 
distributions over different hyperfine sublevels and inhomogeneous spatial distributions. Third, because of low 
mean velocities of atoms, interference effects can be observed during multiple scattering despite the absence of 
any ordering in the atomic cloud. 

Depending on the concentration n of atoms in a cloud, the cloud can interact with light either in the weak 
or strong localization regimes. The weak localization regime assumes that the medium is rarefied so that both 
the mean interatomic distances and the photon mean free path greatly exceed the light wavelength A In this 
case, individual atoms scatter radiation independently and multiple scattering can be treated as a number of 
\^ , successive scatterings by free atoms. By now, this regime has been studied in detail, both experimentally and 



.t^; theoretically [13]- [III. 



The interaction of light with dense atomic clouds has been studied far less thoroughly than with rarefied 
clouds. This is mainly caused by considerable technical difficulties associated with producing such ensembles 
under laboratory conditions. Some problems also arise in the consistent quantum mechanical description of 
their interaction with the electromagnetic field. Under conditions when the mean distance between atoms 
becomes comparable with the light wavelength, it is necessary to consider the resonance interatomic dipole- 
dipole interaction. Because of this interaction, individual atoms cannot be treated as independent sources of 
secondary waves. Conditions are produced for the so-called dependent or cooperative scattering, and we are 
in fact dealing with the interaction of light with a giant quasimolecule combining all atoms in the ensemble. 
Despite these difficulties, the strong localization regime has been attracting considerable attention recently 
because there is good reason to believe that interesting quantum processes are possible in this regime, such 
as the Anderson (strong) localization of light [l^-[l3| and lasing in a disordered medium [23|-[23], which can 
have practical applications. The aim of recent studies in this field was to obtain these effects experimentally 
and describe them theoretically. 

The present paper is devoted to the theoretical study of light scattering by dense atomic clouds under 
conditions of strong localization, i.e., to developing a theory of dependent scattering. Obviously the problem of 
a dense ensemble belongs to the field of macroscopic electrodynamics. It is known that such problems can be 
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solved by using either macroscopic or microscopic approaches. In the former case, description is performed by 
using such averaged characteristics as the field strength and atomic polarization, while the microscopic approach 
is based on the concept of the discrete structure of matter consisting of individual atoms or molecules. The 
macroscopic approach, which we employed to solve the scattering problem, was described in detail in [25^. In 
the present paper, we develop the consistent quantum-mechanical description of the strong localization regime. 
We consider the scattering of quasiresonance radiation by an ultracold disordered ensemble of atoms under 
conditions when the mean number of atoms = in the volume can become comparable with unity. We 
assume that probe radiation is weak and restrict ourselves to the case when the field subsystem contains no 
more than one photon. Atoms are assumed immobile. 

A key problem in describing scattering in this case is analysis of the evolution of an atomic ensemble after 
its excitation by initial radiation and investigation of possible decay channels for this radiation accompanied 
by the creation of secondary scattered photons. In this connection, the given problem is closely related to that 
of collective spontaneous decay in a dense polyatomic system. Moreover, this problem can be reduced to that 
of spontaneous decay by considering the evolution of an atomic subsystem in which only one atom located far 
from the cloud is excited in the initial state. By calculating the dynamics of such a system, we automatically 
describe how a photon emitted by this individual atom is scattered in the cloud. 

The consideration of the problem in this formulation allows us to use methods developed for solving problems 
of superradiance of a singly excited ensemble. How dipole-dipole interaction influences collective spontaneous 
decay in a finite system was first studied for the case of two atoms in (26j and (27j . A consistent quantum- 
mechanical approach developed in these papers was later generalized in a number of papers to polyatomic 
ensembles (see [2^-[33| and references therein). However, almost all papers on polyatomic systems used only a 
two-level atom model. This prevents the correct consideration of the polarization properties of light and, as a 
consequence, adequate description of the resonance dipole-dipole interaction at small distances. In addition, a 
real disordered atomic ensemble was replaced in calculations, as a rule, by a model of an averaged continuous 
medium. In this case, the role of random inhomogeneities was neglected. In our opinion, these simplifications 
considerably restrict the application of these theories for describing experiments with ultracold gases; therefore, 
we do not use them in our paper. Another important difference of our approach is that it can be used to analyze 
the scattering of coherent pulsed radiation, which is also important for interpretation of experiments. 

Our paper contains four sections, conclusions, and an appendix. In Section 2, we solve the nonstationary 
Schrodinger equation to calculate the amplitudes of a number of states of a combined system including atoms 
and electromagnetic radiation. In Section 3, the calculated amplitudes are used to determine the intensity, 
spectrum, and polarization properties of scattered radiation. Section 4 contains the description and results 
of numerical calculations. In the Conclusions, the main results are briefly formulated and the possibilities for 
further investigations in this field are discussed. 



II. STATE AMPLITUDES 



The state of a closed system consisting of atoms and electromagnetic field can be described by the wavefunc- 
tion ■(/;. This function in the Schrodinger representation satisfies the equation 

ih^ = Hi;, (H.l) 



H = Ho + V = Hf + Y,Ha + V. (II.2) 

a 

Here Ha, Hf and V are Hamiltonians of free atoms not interacting with each other and the free field, and 
their interaction operator, respectively. In the dipole approximation used here, we have 

y = -^d('^)E(r,), (II.3) 

a 

where E(r) - is the electric displacement vector operator [4l| . 



E(r) = E(+) (r) + E(-) (r) = z ^ ^ — ^ekaflka exp(ikr) + h.c. (II.4) 

E^^^— are the operators of positive and negative frequency components; ata - is the photon annihilation 
operator in a mode with wave vector k and polarization a; V is the quantization volume; is the dipole moment 
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operator of the atom a, e^a are polarization unit vectors satisfying the transversabihty condition; 

a 

We will solve equation pi.l|) by using the general approach developed in Heitler's papers and considered in 
detail in [s^. We will seek the wavefunction -ip in the form of the expansion in a set of eigenstates {|/)} of the 
operator Hq: 

V-^^foiWIO- (11.6) 
I 

Here, the subscript / defines the state of all atoms and the field. 
Equation (jll.ip gives the system of coupled equations 

3 

for the state amplitudes bj{t). Here Ei is the energy of the lowest state of the system of free atoms and field. 

System (IH.7I) contains an infinite number of equations and cannot be solved in the general case. In this 
paper, as mentioned above, we will study the evolution of only weakly excited states of the combined system. 
In addition, we will calculate all the observables of interest, namely, scattering cross sections, decay constants, 
levels shifts, etc., correct to the second order in the fine structure constant. The higher-order radiative correc- 
tions will not be calculated. In this case, the number of \l) states taken into account in calculations can be 
restricted. As shown in p6| . in this approximation it is sufficient to consider only the states 

V'g = \9,9,--9) f^l^a); Eg ^ hujk; 

V'eo = \9,9,--9,e,g,...g)(E)\vac); Ee^^hoja- (11.8) 
i'e^e, = \9,...9,e,g...g,e,g,...g)(g)\ka); Ee^et=H^a+^^b) + ^vJk■ 
^leTe, uja is the transition frequency of the atom a and the ground state energy of the atoms is set equal to 
zero. 

Note that, by choosing the states considered in the calculations, we go beyond the rotating wave approx- 
imation because we included states with energy h{uja -|- Wf,) -I- hujk, which contain, along with a photon, two 
excited atoms a and b. The consideration of nonresonance states proves necessary to describe correctly the 
dipole-dipole interaction of atoms separated by small distances less than the light wavelength. Indeed, the 
dipole-dipole interaction can be described as the result of exchange by resonance photons. Such an exchange 
can occur in two ways, either via resonant or nonresonant intermediate states. In the first case, an excited 
atom emits a photon, which then is absorbed by an unexcited atom. In the second case, at first, an atom is 
excited and emits a photon and then the initially excited atom undergoes a transition to the ground state and 
absorbs a photon. According to the energy-time uncertainty relation, during the propagation time oi t — r/c 
a photon between two atoms, virtual states can exist with an energy that differs from the resonance (initial) 
energy by the value AE ~ fic/r. When the interatomic distance r is large, it turns out that the second way 
is strongly nonresonant and therefore is usually neglected. If the interatomic distance is smaller or close to 
the reciprocal wave vector r < A ,then t < X/c — 1/w and AE > hck = hut. Therefore, we should take into 
account the second energy exchange channel as well, i.e. the channel involving the intermediate state ipe^cf 
The applicability of the rotating wave approximation in problems of the spontaneous decay of excited states of 
polyatomic systems is considered in more detail, for example, in (ssj -(s^. 

Below, we will consider, along with the spontaneous decay and scattering of one-photon radiation, the 
interaction of an atomic ensemble with weak coherent radiation. It is known that weak coherent radiation can 
be quite accurately represented as the superposition of one-photon and vacuum states. Therefore, we will add 
the state 

i^g' ^\9,g,-g,g)(^\vac); Eg,=0. (II.9) 

to group (jll.Sp . Further simplifications will be performed for a system consisting of -I- 1 atoms. We assume 
that N atoms forming a cloud are identical and can be described by a four-level system with one ground and 
three degenerate excited Zeeman sublevels (the J = J = 1 transition). One of the four-level atoms, 
which we will call a source, can differ in the transition frequency and the excited-state lifetime. The general 
expressions obtained below will be also valid in a particular case when all A^ -I- 1 atoms are identical and form 
one cloud. 

Consider the dynamics of a combined atom + field system, assuming that all atoms at the initial moment, 
except the source, are in the ground state, while the field is in the vacuum state. The source atom is in the 
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superposition of the ground and excited states. This means that at the initial moment two states are excited, 
one of which is g' , and the second state is denoted by o. If bo is the initial amplitude of the o state and bo' is 
that 

of the g' state (|6o|^ + l^o' |^ = 1), then system of equations (|II.7p for amplitudes should be supplemented by 
the initial conditions 

6,(0) - 0; 6o(+0) = bo] V(+0) = K'] i + o,g' . (11.10) 

The solution obtained in this case describes the state of the atomic- field system for positive times i > 0. 

Obviously the amplitude of the state t\}gi (III. 91) in the approximations used does not change during the 
evolution of the system. This is explained by the fact that transitions to this state from other states taken into 
account in the calculations are impossible. The transition from g' to any of the states pi.Sp is also impossible 
(virtual transitions to a state with one photon and one excited atom, which cause Lamb shifts, are taken into 
account by the renormalization of the transition frequency; see the Appendix) , and therefore bgi (t) = bo> ■ 

To determine the time dependence of other amplitudes, we extend, by the Heitler method, the solution of 
system pi.7[) to the negative time semiaxis, assuming that 



bi{t) = 0; i < 0. 

The function bo{t) has a jump for t — 0. We will take this into account by adding the inhomogeneous term into 
the equation for bo{t), 

ih^^-Eibi[t) = ^l^,6,(t) +zMj„(5(t)6„. (11.11) 

3 

where 5ik is the Kronecker delta. The system obtained in this way is solved by passing to the Fourier transform 

bi{t) = j -^bi{u})cxp{-iu}t). 

— oo 

In this case, 

{hu - El) biiuj) = ^y^j(^) + ^^^lobo- (11.12) 

j 

We will separate from (jll.121) the equation for bo, taking into account that Voo = 0. 

(hiu ~ Eo)bo{uj) ^ ^Ki&j(t^) +i^6o. (11.13) 

{huj-Ei)bi{uj) = Viobo{uj)+J2Vi3b3{uj). (11.14) 

We introduce new variables ui{uj) 

bi{u;) ^ ui{uj)bo{uj)^ (fuj - El) . (11.15) 

By introducing a singular function <7 (x) = lim (1 — exp(ikx)) / x , it is possible to prove in the general case [sB] 

that the obtained solution will satisfy the initial conditions. 
Functions ui{lo) satisfy the system of equations 

Ug{uj) = Vgo + VgeUe{uj)(. {hu> ~ Ee) , 

Uee(^) = Vee;o -^^^Vee;eUe{uj)(^ {hw ~ Ee) , (11.16) 

67^0 

Ueiuj) = YVegUg{uj)<, (huj - Eg) +YVe-eeUeei^)'^ {fl^ - Eee) , e ^ O, 
g ee 

By solving this system, we obtain the amplitude bo{uj) 

bo{u;) = ^VtT' (II-l^) 

So(w) = A - tr/2 = ^ VogUg{uj)c: {hcj - Eg) /h+Y, Vo,eeUee{uj)'; {tlLJ - E,e) /fl (11.18) 
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The subscript of the excited atom is omitted for brevity in relations (jII.16P - (jII.18[) . and it is also taken into 
account that the relation (x) = 1 is fulfilled for singular functions ^ (x) . 

By using the first and second equations of system (IIL16p . we obtain the closed system of equations 



.{uj) = ^ VegVgo<, (hhJ - Eg)+'^ Ve-eeVee-o^ {fVjJ - E'ee) ' 



VegVge'^ {fiLO - Eg) + ^ Ve-eeVee;e"i {fl^ " -^ee) 



Ue'{uj)<;ihuj - Ee') , o. (11.19) 



for the states with one excited atom in the cloud. 

Note that this system, unlike (jll.lll) . contains a finite number of equations. This number is determined by 
the number of atoms in the cloud and the number of excited states considered in the calculations. The sums 
X] ^egVge'<; {fujj — Eg) and ^ Ve;eeV'ee;e"^ — E^e) Can be Calculated for an infinitely broadband field reservoir 

g ee 

by assuming that the quantization volume is infinite (see the Appendix). By using expressions presented in 
the Appendix, we rewrite system pi. 191) in a more compact form: 

[(W - iOa)5ee' " ^ee'{l^)] Ue'{uj)<i {flU ~ E^') = 5]eo(w) (11.20) 

Here, the vector I]eo(a;) in the right-hand side determines the excitation of different atoms in the cloud by the 
source radiation, while the matrix Y^ee' (w) describes radiation transferred between atoms in the cloud. 

System pi. 201) can be reduced to an integral equation by using the continuous medium approximation. This 
significantly simplifies the solution of the problem for a two-level atom system [s^-js^]. Moreover, in this 
case, even an analytic solution is possible for spatially homogeneous spherical clouds. This solution neglects, 
however, the important properties of real physical systems, and therefore we will solve the linear system pi.20[) 
numerically. The numerical solution is restricted by the possibilities of available computers which allow us 
to consider clouds containing a few thousands of atoms. The number of atoms in quasi static traps, where 
the required concentrations can be achieved, can be two orders of magnitude greater; however, by using this 
approach, we can correctly describe all polarization effects taking into account random inhomogeneities of the 
medium. Both these factors, as follows from our previous calculations 25|, are quite important. 

The formal solution of system pi.20p can be written in the form 

Ue(w)'? {flUJ ~Ee)^Yl Ree'{uj)Ye'o{i^), (11.21) 

where 

R,,,{lj) = [iu~LUa)S,e' " See'(w)]"' ■ (11.22) 

By using pi.21l) and pi. 221 we obtain the amplitudes of all states of interest for us: 

bg,{t) = b,r, (11.23) 

OO 

Kit) ^ ^ ^oexp(-^^) 

— OO 



m = I ^ ^^Vn ' ("-2^) 

J 271 UJ - UJo -2^o[^) 

— OO 
OO 

5ee(t) = / l!^ K^i^~E,^)-M~^-^) I y^^,^ + Y: V.e.M^)^M^) 1 ; (11.27) 
/ ZTT id — UJo ~ ^o(^) * ' 



-OO 



T.o{uj) = Soo(a;)+ ^ I]oe(a;)i?ee' (w)Se'o(w) (11.28) 

Relations pi.23P - pi.28[) completely describe the studied system under the above assumptions. In the next 
section, we will use them to obtain explicit expressions for scattered radiation parameters observed in exper- 
iments. However, an important preliminary remark should be made. Expressions pi.23p - pi.28|) are valid for 
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any distance between an atom excited from the ground state and the cloud. They describe, in particular, 
the modification of the spontaneous decay of the atom near the cloud or even inside it. Our main goal is to 
describe the scattering of light by a dense atomic cloud. Therefore, below we restrict ourselves to the case when 
a source atom is located far from the cloud. In this case, we can neglect the reverse influence of the cloud on 
the spontaneous decay of this isolated atom taking into account only the first term in (jll.28p by calculating 
So(w). By including the Lamb shift in the definition of frequency Wo, we obtain 

Eo(w) = -i7o/2. (11.29) 

Here, 70 is a decay constant of the atom that was initially excited. This constant can be arbitrary. 

Along with simplification of the quantity 'E,o{uj), the expression for the vector Eeo(ti;) describing excitation 
of different atoms in the cloud by the source is also simplified in this case. Because the distance from the cloud 
to the emitting atom is large, retardation effects are important and the polar approximation cannot be used. 
However, we can use the rotating wave approximation and retain only one term decreasing the least rapidly 
with increasing distance from the source to the cloud. In this case 

- - E -i-^ ~ [-) ^-p j ; (11-30) 

l^LM L.U J 

Here, r^o is the vector drawn from the source oo to the atom excited to the e state (hereafter for brevity, we 
will call it the e atom). We will transform expression (jll.30p taking into account that Tgo greatly exceeds the 
size of the system. We introduce a coordinate system with a center inside the cloud. Let re an be the 

radius vectors of the cloud atoms and the source, respectively. Then, roe = \J (fo — fe)^ — 'i'o{l — ToVe/rl) = 
To - nr^; n = r^/ro and 



Eeo(w) -E 



hro 



exp (ifcro — ikre) ; (11.31) 



Here, k =a;n/c. We can see from ()II.31|) that the source radiation in the atomic ensemble region can be treated 
quite accurately as a plane wave. 

III. TRANSITION PROBABILITIES AND SCATTERING CROSS SECTION 

Knowing the amplitudes of the states, we can calculate the probability of finding the system in a specified 
state at any instant of time. In experiments, as a rule, the parameters of electromagnetic radiation are 
measured at large distances from a scattering cloud. These parameters are determined by the quantity bg(t),. 
The calculation of this quantity is the main purpose of our paper because different states g correspond to 
different frequencies, polarizations, and wave vectors of the emitted photon. Relation (|II.25p allows us to find 
bg{t) at an arbitrary instant of time. Below we consider the two most important particular cases allowing a 
significant simplification in the analytic form. The first case is the calculation of the probability of different 
field states after a long time interval when spontaneous decay has ended. In our case, these are times exceeding 
the excited-state decay times of the source atom and atomic ensemble. The second case is the calculation of the 
transition probability per unit time when the system state still only weakly differs from the initial state. These 
times are considerably shorter than the decay time (but, of course, they greatly exceed reciprocal transition 
frequencies) . 

We will solve this problem taking lim [hw — Eg) exp(— iw<) = —2niS{huj — Eg) exp(— ii?„t//i). 



^ait) = 1 /'°"^P( '^^1^ ^ \Vgo+ J2 VgeRee'{u^g)^e'o{u;g) ] , t ^ +00 (IH.l) 



Here ujg = Eg/h is the frequency of the emitted photon. 

The first term in expression (jlll.ip describes radiation incident directly on a photodetector from the initially 
excited atom. The angular and spectral distributions of radiation from the cloud are determined by the square 
of the modulus of the second term in parentheses in (IIII.ll) . We will assume here that the photodetector is 
equipped with a set of apertures eliminating the contribution of the initially excited atom to the photocurrent. 
Then, the probability of finding a photon in the final state in a given mode will be 



P, 



bo J2 VgeRee'{i^g)^e'o{uJg) 



' {Cg^^of+l'o/^ 



(111.2) 
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The dominator in this expression is related to the efBciency of excitation of the cloud by the source radiation. 
We calculated it using expression pL29[) . For small 70, only a narrow spectral group of the cloud states is 
efficiently excited. To determine the total excitation spectrum of the cloud in this case, it is necessary to repeat 
experiments many times by changing the source radiation frequency Wo- We can obtain the cloud radiation 
spectrum in one experiment by using a radiation source emitting a very short pulse, so that 70, is large and all 
the characteristic collective excitation frequencies lie within the pulse spectral width 70,. 

In another limiting case, when 70, is much smaller than the characteristic rates of transient processes in the 
cloud, a quasi-stationary process will take place during which the quasi-stationary distribution of excitation 
over the cloud will be established. This distribution will then change slowly (with the decay rate 7q"^ of 
the excited state of the source atom). In this case, we can calculate the transition probability per unit time 
w = d \bg{t)\ jdt. To do this, it is necessary to perform the passage to the limit 70 before passing to large 
times: 



70-i-O 



iboduj 

2TTh 



^ (w - Wg) <J (w - Wo) exp(-iwi) I + ^ VgeRee' ii^)^e'o{^) \ (III. 3) 



We will calculate w using the expansion of the product (w — ujg) <; (w — Wo) into simple fractions: 
q-yg {uj — cjo) W (w — Wg) = [<^7o ("^ ~ ^o) — 'fcr (w — ujg)] (^o — ujg) . Here a and 70 are constants. By having 

them tend to zero, we obtain the corresponding generalized function (further we will assume for definiteness 
that cr I 7o, which does not affect the final result). Taking into account that (wo — Wg) — <r*_^^ (wo — Wg) = 

—2TTiS {ujo ~ LOg), we obtain after simple transformations 



w — 2t: \bo\ 



ygo+ E VgeR ee' ('^o)^e'o('^o) 



-5 [uJo - UJg) ■ 



(III.4) 



As earlier, the first term describes the contribution of radiation emitted directly by the initially excited atom. 

Expressions (|nL2|) - (|nL4t can be used to describe experiments in which the occupation of different modes 
of the field is detected, for example, when modes with a narrow group of wave vector directions are selected 
with the help of a convergent lens. To describe experiments without such a selection, for example, with point 
photodetectors, it is necessary to know the mean parameters of the field at a certain spatial point at a specified 
instant of time. These mean parameters can be found from the known wave function 



(III.5) 



For example, the mean projection a of the vector E'^"'^^ at the point r at the instant t is 



3:E(-)(r) 1^) = / ^:^ ^bob:,eM-^.t) U ^ £^^(.)^^,(^)s,„(.) ] . (III.6) 



Here, the matrix T,aei^) describes the propagation of radiation from an atom excited in the state e to the 
observation point. 



E„e(^) - J2 ("'^ <E(+)(r) \g) Vge^ {hw - Eg) /h. 



(111.7) 



The explicit expression for this matrix can be obtained as it was done in the Appendix upon calculation of 
infinite sums entering equation ()II.19|) . In the rotating wave approximation, we obtain 



1 



h\r - r,. 



« (r - re)) (dg.e (r - r^)) 



r - re 



exp I 



w r — r„ 



(HI. 



Here, re is the radius vector of the atom excited to state e and r is the radius vector of the observation point 
in the chosen coordinate system. 

The polarization components of the intensity correlation function can be calculated similarly: 



,E(^)(r)e:E(+)(r)|^) = |feof 



00 / 

f iduj hexp(-iujt) I ~ ~ , , , , 

/ 1^ r~- 77; \^ao[UJ)+ y ^ Sae(w)-ftee'('^)^e'o('^) 

J 2-K U -UJo + l^o 2 \ ^ 



(IH.9) 
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Expressions pil.6p and pil.9[) are considerably simplified in the quasistatic case. By retaining only the 
contribution of light scattered by the cloud and passing to the limit 70 — 

(^|e*E(+)(r)|^) = hbX'eM-i^ot) ^ E„eK)i?ee' K)Se'oK); (III.IO) 

e,e' 

2 

(^|e„E(-)(r)e:E(+)(r)|V;) = {bof S„eK)i?ee' (wo)Se'oK) • (ULll) 

Knowing the transition probability per unit time, we can determine the differential cross section of light 
scattering by an atomic cloud. We calculate this cross section by assuming that the distance from the source to 
the cloud is so large that a wave coming from the source to the atomic ensemble is a plane wave with the wave 
vector k =ujon/c with good accuracy [see (jll.3ip ]. This wave is polarized. The direction of its polarization unit 
vector e is determined by the orientation of the vector k and a Zeeman sublevel to which the source atoms was 
excited. We will find the photon flux density J of this wave in the atomic ensemble region from relation (jIII.9|) 

by retaining only the first term in it: J = \bo\ h T,ao{^o) /2ttuJo- Here, the summation is over any two 

a 

basis polarizations in a coordinate system determined by vector k. 

To obtain the explicit expression for the scattering cross section, we retain in expression (IIII.4[) for the 
probability w only the contribution caused by the cloud, sum it over a narrow group of finite states of the 
field, and divide by the fiux density J of photons incident on the atomic ensemble from the source. Taking into 
account that the radiation frequency does not change upon scattering in this case, we obtain the expression 
k',e' 



da 
dJ7 



^ {e'*dg.e) Ree' (Wo) {^de'-g) CXp {i{kre' - k'Ve) 



(III.12) 



for the cross section of photon scattering from the mode k, e , to an arbitrary mode k', e'. 
Note that we obtained the relation for the total cross section earlier |23] by using the T-matrix formalism 
and showed, in particular, that the Ree' matrix is the projection of the resolvent of the system under study 
on the states with one excited atom. The relations obtained in this section will be further used to analyze the 
spectral, angular, and polarization characteristics of light scattering by ultracold atomic clouds. 



IV. RESULTS OF CALCULATION 



The main quantity necessary for the description of scattering is a resolvent Ree' (w) because it determines 
amplitudes (|II.23|) - (|II.28I) of all quantum states involved in the problem and all experimental characteristics of 
the process under study [see, for example, (jIII.6|) . pil.9[) and pil.121) '). In our paper, this resolvent was found 
numerically. Depending on convenience, we calculated either the reciprocal matrix pi.22[) or directly found the 
numerical solution to system pi.20p . The resolvent value considerably depends on the spatial configuration 
of the ensemble under study. The analysis of random inhomogeneities of a medium and averaging over these 
inhomogeneities were performed by the Monte Carlo method of statistical tests. The spatial distribution of 
atoms in the cloud was specified by means of pseudorandom-number generators, and then a particular observable 
of interest for us was calculated. This observable was calculated repeatedly for different random configurations, 
and the results were averaged. Note that in most of the papers studying collective effects in dense media, the 
averaging procedure is performed differently by replacing the real medium by the averaged medium with a 
continuous density distribution, and then the observables are calculated for this averaged medium. As follows 
from our analysis [25i] , in this case, some correlation effects affecting observables are described incorrectly. 

The expressions obtained in the previous section describe light scattering by ensembles of ultracold atoms 
under various experimental conditions, in particular, for spatially inhomogeneous clouds of different shapes 
and arbitrary densities. We are mainly interested in dense clouds with high atomic concentrations in which 
dependent scattering effects are important. However, we will not restrict ourselves to this case only. Some 
results obtained in the paper concern optically dense but rarefied atomic ensembles. The analysis of the weak 
localization regime as a limiting case is important, first, for verifying and substantiating the method proposed 
and, second, it allows us to refine some results that were earlier obtained in this regime by other methods. 

Below, we analyze the two main characteristics of scattering measured in experiments, namely, the angular 
(spatial) distribution of scattered light and the spectral dependence of the scattering cross section. We restrict 
ourselves to the results obtained for the quasistatic case. 
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A. Spatial Distribution of Scattered Radiation: The Effect of Coherent Backscattering 

The scattering of light by an opticahy dense medium in the presence of locahzation phenomena is considerably 
determined by the effects of interference and radiation trapping. Therefore, we verify first of all that the 
approach proposed in this paper correctly describes these effects. 

The basic equations (jll.20[) do not contain derivatives with respect to spatial variables. In this connection it is 
interesting to study how these equations describe radiation transfer in an optically dense medium, in particular, 
how they reproduce the Bouguer-Lambert law. According to this law, the intensity of the coherent component 
of radiation transmitted through a medium layer is attenuated exponentially. The exponent b — naL (the 
optical thickness of the layer) for a homogeneous medium is determined by the concentration n of atoms, the 
layer thickness L, and the scattering cross section a by an atom. For the J = •H- J=l transition in a rarefied 
medium, we have a = dirX^. 

To analyze the radiation transfer process, we consider the scattering of resonance radiation by a model 
cylindrical atomic cloud with radius R and length L. The wave vector of the incident wave is assumed parallel 
to the cylinder axis. We will calculate the spatial intensity distribution Icoh for the coherent component [which 
is determined by the square of the amplitude modulus (jIII.6[) )] in a plane perpendicular to this axis. The 
distance from the rear of the cylinder to the observation plane is chosen equal to a few A to fit several 
Fresnel zones within the end for any point of this plane. The results of calculations for a rarefied medium with 
nX^ = 10~^, R = 24A, Zd = 4A, and different L are shown in[T] 

The obtained intensity distribution corresponds to the Fresnel pattern for diffraction by a semitransparent 
screen. Outside the shadow region, oscillations are observed typical of diffraction from the screen edge. The 
intensity attenuation in the geometrical shadow region is caused by the destructive interference of the source 
radiation and secondary radiation emitted by the atomic ensemble. This interference formally appears in our 
approach when two terms in (|III.6|) are taken into account. Note that the diffraction pattern in our approach 
is obtained directly upon summation of contributions from a great number of microscopic scatterers forming 
the medium. 

The transmission coefficient Tcoh was calculated by finding the mean light intensity in the geometrical shadow 
region. These mean intensities for two thicknesses are shown in[l]by the dashed lines. The dependence of the 
sample transmission on its length L showed that the Bouguer-Lambert law is reproduced quite accurately in 
our calculations. In this case, the optical thickness increases linearly with increasing L. As for the dependence 
of optical thickness b on concentration n, it increases linearly only up to concentration nA"^ < 10^^ ([2]). The 
deviation from linearity at higher n is caused by collective effects, namely, by the interatomic resonance dipole- 
dipole interaction. This interaction leads to shifts of atomic levels, so that the scattering cross section a for 
resonance radiation decreases. This decrease partially compensates the increase in concentration: the higher 
the atomic density, the better the compensation. This effect, in particular, can partially explain the difficulty 
in experimental observation of strong localization in atomic gases. 
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Figure 1: Spatial intensity distribution Icoh of the coherent component of light scattered by a rarefied cylindrical 
cloud. Coordinate x is measured in a plane perpendicular to the cylinder axis. The cloud radius is _R = 24A and the 
concentration is nX^ — 10~^. The different curves correspond to different optical thicknesses 6 = naL of the scattering 
cloud. 
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Figure 2: Transmission coefficient of the atomic cioud for the coherent component as a function of cloud concentration 
(solid curve). The deviation from the Bouguer-Lambert law a = 6nX^ (dashed straight line) on passing from rarefied 
to dense clouds is clearly seen. 

We now consider light scattering in the case of Fraunhofer diffraction. We analyze the angular distribution 
of the cloud radiation at a large distance from it. The calculation is performed using expression (jIII.12[) for the 
differential scattering cross section. We assume that the light frequency coincides with the transition frequency 
in a free atom, while the probe radiation is circularly polarized. We will analyze, in particular, coherent 
backscattering (CBS); therefore, to avoid problems with sharp boundaries, the atomic cloud will be assumed 
inhomogeneous with a Gaussian spatial density distribution. 

Figure [3] shows the angular distribution of scattered light for two polarization channels, with light helicity 
preserved (H || H) or changing to the opposite (H ± H). The calculations were performed for a cloud with 
a Gaussian radius R = 12A.. The density at the cloud center was nX^ = 10~^. This corresponds approx- 
imately to the intermediate value between dense and rarefied ensembles. When the helicity was preserved, 
we observed typical Fraunhofer diffraction. The principal diffraction maximum caused by coherent forward 
scattering dominates. Several higher-order peaks are well distinguished, although because of a smooth decrease 
in the concentration at the edges of the Gaussian cloud, these peaks are not as distinct as for diffraction from 
objects with sharp boundaries. Along with diffraction peaks, a comparatively weak background is observed, 
which is caused by multiple scattering in the cloud. The background intensity weakly depends on angles, except 
for a narrow region near backward scattering, where a local maximum is observed. This is the CBS cone. 

Radiation with orthogonal helicity (the dashed curve in Fig. |3]) appears due to incoherent multiple scattering 
(incoherent light scattering by different atoms). The intensity of this radiation in the case under study is low. 
It increases monotonically for scattering angles above tt/2. The CBS effect for this polarization channel and 
the transition under study is known to be a very weak (see, for example, jTEj)) and is almost unnoticeable 
at this scale, which is correctly reproduced by our calculations. Thus, our microscopic calculation reproduces 
the interference effect upon multiple scattering in a disordered medium for both main scattering channels. 
However, unlike methods used earlier for calculations of this effect and taking into account only a finite number 
of scattering orders, our approach takes into account all processes of multiple scattering of photons in a medium 
and also allows us to determine the cone shape at arbitrary atomic concentrations. 

The CBS properties are shown in more detail in Fig. |4l This figure demonstrates the concentration depen- 
dence of the CBS cone for the H || H channel. The visible shape of the cone changes comparatively weakly in 
this range of n. Figure 4 demonstrates the qualitative change in the CBS effect for dense clouds. It is known 
that all weak-localization theories proposed so far predict the maximum value of the amplification factor (the 
ratio of the cone maximum to the value of the quasi-isotropic pedestal) equal to two. This is explained by the 
fact that the cone amplitude is caused by the interference of the waves scattered by the same chain of atoms 
but propagating along this chain in opposite directions. These waves in rarefied media interfere only in pairs. 
We see that our calculation for rarefied media (the lower curve in Fig. IH) quite accurately gives an amplification 
factor of 2. The amplification of backscattering in denser media exceeds 2. In our opinion, this is explained by 
two main reasons. First, in dense media, along with the conventional mechanism resulting in weak localization, 
the interference of beams scattered by different but spatially close chains makes a certain contribution. The 
second reason is the decrease in the fraction of light scattered outside the cone at comparatively large angles, 
i.e., the decrease in the background intensity with respect to which the amplification effect is measured (see 
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Figure 3: Angular distribution of scattered light intensity for two polarization channels. The solid and dashed curves 
correspond to scattering without a change in the helicity (H || H) and with a change in helicity to the opposite (H _L 
H). 
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Figure 4: Cone shape for CBS by a spherically symmetric Gaussian cloud with a radius of 12X as a function of the 
atomic concentration 



B. Scattering Spectrum 

We consider the spectral dependences of the differential and total scattering cross sections as another example 
of the application of the theory developed above. The total cross section can be obtained by integrating the 
differential cross section over the entire spherical angle or, according to the optical theorem, by calculating the 
amplitude of light scattered at the zero angle. In this paper, we used both these methods, the optical theorem 
being fulfilled independently of the cloud density. 

Figure [5] shows the total scattering cross section spectra for a spherically homogeneous, on average, cloud with 
radius R = 15A at different densities. At low concentrations, when the interatomic interaction is insignificant, 
we see a monotonic spectral dependence, which is typical of scattering by an isolated atom. The spectrum is 
somewhat broadened compared to the monatomic scattering spectrum because this medium is not optically 
thin even at these concen trations. As the cloud density is increased, the shape of the spectrum begins to be 
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distorted. At the concentration nX^ = 0.2, the spectral dependence is substantially nonmonotonic, and several 
distinct local maxima appear. 




Figure 5: Spectra of the total cross section of light scattered by a spherical cloud with radius R — 15X for different 
atomic concentrations. 

This behavior of the cross section can be explained either in the context of macroscopic electrodynamics or by 
using the microscopic approach. Macroscopic electrodynamics describes scattering as the result of multibeam 
interference (or diffraction) from a sphere with the specified permittivity. Such scattering in the case of spatially 
homogeneous spheres with arbitrary permittivities was considered in detail in the Debye- Mie theory (see, for 
example, jssj). This theory showed, in particular, that the maximum scattering cross section for transparent 
spheres with the refractive index n was observed for 2i?fc(n — 1) ~ 4.. In this case, the scattering cross section 
was approximately four times greater than the sphere cross section area ttR^. As the radius is decreased 
compared to R, the ratio of the scattering cross section a to ttR^ tends to zero. As R is increased, the ration 
a/irR^ oscillates and tends to 2 for i? ^> oo. The scattering cross section behaves similarly when the sphere 
material has a lower optical density than the environment. However, there exists some quantitative difference. 
The maximum scattering cross section is lower and is observed for smaller values of R. If the medium is not 
transparent, i.e., the refractive index has a noticeable imaginary part (the absorption coefhcient), the diffraction 
type changes. The oscillations mentioned above are considerably smoothed, although the maximum scattering 
cross section can exceed 2ttR^. 

This behavior of the scattering cross section agrees qualitatively with results shown in Fig. [5] In the case 
of an atomic cloud, the permittivity is complex. We can introduce the real refractive index and absorption 
coefficient (although note that in the case under study, a wave is not really absorbed by an atomic ensemble 
because the escape of photons from the incident wave-mode is accompanied by the appearance of photons in 
other modes). The microscopic approach developed in our paper gives accurate values for the permittivity 
of a homogeneous atomic medium with the specified density. However, the discussion of this issue is beyond 
the scope of this paper. Note that the permittivity value is determined by studying the spatial distribution 
of the atomic polarization taking into account that for detunings A exceeding the atomic resonance width, 
the refractive index decreases with increasing A considerably slower than the absorption coefficient. For the 
specified radius of the cloud Rk = 15, the condition 2Rk{n — 1) ~ n — 1 ~ 0.1333. 

Our calculation of the permittivity shows that the position of the left maximum qualitatively agrees with 
the predictions of the Mie theory. The intensity of this maximum is smaller than the maximum possible value 
for transparent media, which is related to a small but finite absorption coefficient in this frequency region. 
As the atomic concentration in the cloud is increased, the maximum shifts to the region of greater detuning 
magnitudes and its amplitude increases. This is explained by the fact that a comparatively small value of n 
in a denser medium is achieved for a large magnitude of detunings A. The relative influence of absorption for 
large detunings also decreases. 

The maximum observed for positive detunings also can be explained by the Debye-Mie theory. The refractive 
index of the atomic medium in this frequency region is smaller than unity, and therefore the right maximum 
intensity is lower and it is located at smaller magnitudes of A than the left one. 

The absorption coefficient for small detunings is large compared to n — 1 (the difference of the real refractive 
index from unity), and therefore there exists another peak related to the absorption coefficient maximum. This 
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peak is most distinctly observed at large concentrations when the spectral positions of all three cross section 
maxima are well distinguished. 

Figure [5] shows the concentration dependence of the shape of the total cross section spectrum for a fixed 
cloud size. Figure [S] illustrates the dependence of the cross section on cloud size for the specified concentration. 
The calculation was performed for nX^ — 0.3 when the interatomic interaction is considerable. We can see that, 
as the cloud size increases, the maxima caused by refraction shift to the region of large detunings, according 
to the condition n — 1 ~ 2/Rk. The position of the central maximum at the specified concentration does not 
change. 
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Figure 6: Spectra of the total cross section of light scattered by spherical clouds of different sizes. The atomic concen- 
tration is nX'^ = 0.1. 

The frequency dependence of the total cross section can be explained by using the microscopic approach. To 
do this, it is necessary to consider the distributions of collective atomic states over energies and lifetimes, more 
exactly, the distribution of the poles of the resolvent Ree' (i^) over the values of their real and imaginary parts 

In light scattering experiments, the differential cross section is most often studied rather than the total cross 
section. As a rule, the intensity of light scattered through a specified angle is measured [l^ and its dependence 
on the probe light frequency is determined. Figure [7] shows this dependence calculated by us for a spherical 
cloud with radius Rk = 15 and concentration nk~^ — 0.1. On the ordinate, the total light intensity (without 
selection over polarizations) is plotted. Different curves correspond to scattering through different angles ?. 
Note that the behavior of these curves is substantially different. For some scattering angles, the spectral 
dependence is monotonic, whereas for others, the nonmonotonic behavior is observed. A similar difference will 
be manifested in the fiuorescence of the cloud excited by a light pulse. In the case of a short pulse whose 
spectrum overlaps the entire excitation spectrum of the atomic ensemble, fiuorescence in some directions will 
decay monotonically, while in other directions, beats will be observed. 

The spectral behavior is different not only for scattering in different directions but also for different polar- 
ization components scattered in the same direction. This is illustrated in Fig. [51 which shows the intensities 
of two orthogonal polarization components as functions of the incident light frequency. The calculation was 
performed for the same conditions as in Fig. [71 The scattering angle is 7r/3. The polarization of incident light, 
as in all other calculations, was circular. The results demonstrate a considerable dispersion of polarization: the 
polarization of light scattered in the given direction can considerably depend on the incident light frequency. 

V. CONCLUSIONS 

We have developed a consistent quantum- mechanical theory of cooperative scattering of weak electromagnetic 
radiation by an optically dense macroscopic atomic ensemble. The microscopic approach proposed in the paper 
can be applied both to rarefied clouds, when the average number of atoms in the volume is small, and to the 
clouds in which this parameter exceeds unity. The differential and total dependent scattering cross sections 
have been calculated for quasi stationary radiation using this approach. The angular, spectral, and polarization 
properties of scattered light were determined. The dependence of these characteristics on cloud size and atomic 
concentration was studied. By summing the contributions from many atomic scatterers, we obtained, along with 
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Figure 7: Intensity of light scattered in different directions as a function of the incident radiation frequency. 
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Figure 8: Spectral dependences of the intensity of light scattered through angle 7r/3 for two polarization channels: solid 
curve and dashed curves are scattering without a change in the helicity H || H and with a change in the helicity to the 
opposite H ± H. 

the known Fresnel and Fraunhofer diffraction phenomena, the contribution from multipole incoherent scattering 
and calculated the interference effect of coherent backscattering. The parameters of the backscattering cone 
were determined. Note that, unlike approximate methods used earlier, in which backscattering was calculated 
taking into account a finite number of scattering events, our approach allows us to calculate these parameters 
exactly. This effect was calculated for the first time for dense clouds where the mean free path is comparable 
with the light wavelength. It was found that the complex nature of interference effects, when aside from the 
traditional mechanism resulting in a weak localization, it is necessary to consider the interference of beams 
scattered by closely located chains, leads to an amplification factor exceeding 2. 

The dependence of the scattering cross section on the incident light frequency has demonstrated a strong 
modification of the spectra with increasing atomic concentration. It was shown that these spectra for dense 
media have a prominent nonmonotonic nature and exhibit several maxima. It was also found that collective 
effects lead to a considerable difference in the spectral dependences of scattering in different directions. The 
investigation of the polarization properties of scattered light revealed a considerable dispersion of polarization, 
which strongly depends on the incident light frequency. 

Our study concerns the case when the motion of atoms can be neglected. It is obvious, however, that this 
motion can play an important role in some experiments. Therefore, we believe that it is interesting to consider 
further the influence of atomic motion on collective effects in dense atomic systems and the interaction of light 
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with such systems. In addition, it is important for practical purposes to generalize the microscopic approach to 
the case when an additional control field is present. In our opinion, the use of high density atomic systems to 
obtain electromagnetically induced transparency and to " stop light" can be quite efficient for solving problems 
in quantum informatics. 
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VI. APPENDIX 



Infinite sums ^ VegVge' <^ {fiuj — Eg) and ^ Ve-eeVee;e"^ {ti^ — E^e) appearing in Eqs. pi.l9p can be calculated 

g ee 

for an infinitely broadband field thermostat assuming that the quantization volume is infinite. We first consider 
the first sum, assuming that in state |e') the atom a is excited, while in |e) atom h is excited. In dipole 
approximation (|II.3p . we have 
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cdk sin(A:r) 
nr uj ~ kc + 



(VI.l) 



Here, r — r;,— r^. Expression (jVI.l[) was derived by calculating the sum over polarizations using completeness 
relation pi.5p . replacing the sum over k by an integral, and using one of the representations of the singular 
function <;{x). 

Formally, some of the integrals in (|VI.1|) diverge. However, this divergence is related to the behavior of 
integrands for large k, which is caused by the use of the dipole approximation. By introducing, as usual, a 
slowly decaying exponential factorexp(— efc) with a small exponent s > 0,, we can provide the divergence of 
integrals. By calculating integrals and tending the parameter e to zero in finite expressions, we obtain 
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(VI.4) 



Here, Ci{x) and Si{x) are integral cosine and sine, respectively. 

Similarly, we can calculate the two-atom contribution to the second sum caused by nonresonancc intermediate 

states 
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(VI.5) 



Here Wq is the transition frequency of atoms forming the cloud. It is easy to see that this contribution for small 
r is comparable with resonance term (jVI.2p . 

Expressions (jVI.2p and (jVI.Sp can be considerably simplified in the so called polar approximation, when their 
value for frequency uj is replaced by their value for frcquencywa of the atomic resonance. This approximation 
was studied in detail in [40[ , where it was shown that it can be applied in systems where retardation effects are 
insignificant. The value of R for atomic clouds in quasistatic traps, which are studied in this paper, does not 
exceed hundreds of micrometers. The shortest evolution times of the atomic system to be compared with the 
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delay times are the lifetimes of short-lived superradiance states. For systems with size R, these times are on the 
order of T{nk~^R)~^ where r is the natural excited state lifetime of a free atom. Under such conditions and 
for typical parameters of alkali metal atoms, retardation effects can be neglected even for clouds with density 
nk~^ ^ 1. In this case, using the polar approximation, we obtain 




(VI.6) 
(VI.7) 



Note that the latter expression can be also rewritten in terms of spherical Hankel functions and spherical 
harmonics [25ll3T|. 

The monatomic contribution, when e' = e, has the form 




(VI.8) 
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where 7 is the spontaneous decay rate, which is the same for all the sublevels of the excited state. 

The sum '^Ve-^eeVee-.e'i {fy^ — Eee) gives Only the Lamb shift of the ground state because it describes the 

ee 

spontaneous excitation of an atom located close to the excited atom, which is accompanied by spontaneous 
decay with absorption of a virtual photon. We neglect the Lamb shifts of the ground and excited states, 
assuming that they are included in the transition frequency oja- 
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